Multi-plateau magnetization curves of one-dimensional Heisenberg ferrimagnets 
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Ground-state magnetization curves of ferrimagnetic Heisenberg chains of alternating spins S and 
s are numerically investigated. Calculating several cases of (S,s), we conclude that the spin-(S, s) 
chain generally exhibits 2s magnetization plateaux even at the most symmetric point. In the double- 
or more-plateau structure, the initial plateau is generated on a classical basis, whereas the higher 
ones are based on a quantum mechanism. 

PACS numbers: 75.10.Jm, 75.40Mg, 75.30.Kz 



I. INTRODUCTION 

Ground-state magnetization curves of low-dimensional 
quantum spin systems have been attracting much recent 
interest due to their nontrivial appearance contrasting 
with classical behaviors. A few years ago there appeared 
an epochal argument |Q in the field. Generalizing the 
Lieb-Schultz-Mattis theorem H|| , Oshikawa, Yamanaka, 
and Affleck proposed that magnetization plateaux of 
quantum spin chains should be quantized as 



?unit ui = integer . 



(1.1) 



where S nn it is the sum of spins over all sites in the unit 
period and m is the magnetization M divided by the 
number of the unit cells. Their argument caused re- 
newed interest Q] in the pioneering calculations |5|-[7j of a 
bond-trimerized spin-i chain and further stimulated ex- 
tensive investigations into quantum magnetization pro- 
cess. Quantized magnetization plateaux were reasonably 
detected for spin-± |§, spin-1 and sph>§ |||l3) 

chains with modulated and/or anisotropic interactions. 
Totsuka |0, Cabra and Grynberg |^5| , and Honecker 
p6[ developed calculations of general polymerized spin 
chains. Numerous authors have been making further 
theoretical explorations into extended systems including 
spin ladders Jl7]-^0| and layered magnets |pl|-p3|| . Exper- 
imental observations |24|,|25| have also followed. 

Mixed-spin chains, which have vigorously been studied 
in recent years pf^-|3r^|, also stimulate us in this context. 
Theoretical investigations into them are all the more 
interesting and important considering an accumulated 
chemical knowledge on ferrimagnetic materials. Kahn 
et al. |37[] succeeded in synthesizing a series of bimetal- 
lic chain compounds such as MM'(pba)(H20)3-2H20 
(pba = 1, 3-propylenebis(oxamato) = C7H 6 N 2 06) 
and MM'(pbaOH)(H 2 0) 3 (pbaOH = 2-hydroxy-l, 3- 
propylenebis(oxamato) = C7H6N2O7), where the alter- 
nating magnetic ions M and M' are flexible variables 
and therefore the low-dimensional ferrimagnetic behav- 
ior could systematically be observed. Caneschi et al. 



P8| synthesized another series of mixed-spin chain com- 
pounds of general formula M(hfac)2NITR, where metal 
ion complexes M(hfac)2 with hfac = hexafluoroacetylace- 
tonate are bridged by nitronyl nitroxide radicals NITR. 
There also exist purely organic molecule-based ferrimag- 
nets p^| , where sufficiently small magnetic anisotropy, 
whether of exchange-coupling type or of single-ion type, 
is advantageous for observations of essential quantum 
mixed-spin phenomena. 

Magnetization curves of Heisenberg ferrimagnetic 
chains were revealed by Kuramoto p3| . His argument 
covered an effect of next-nearest-neighbor interactions 
but the constituent spins were restricted to 1 and i. Al- 
though an alignment of alternating spins S and s (S > s) 
in a field, as described by the Hamiltonian 

N 

n = Y, [ (1 + s )(Sj ■ sj)* + (! - • S i+ l) a 

i=i 

- H(S* + , (1.2) 

with (S ■ s) a — S x s x + S v s v + aS z s z , is so interesting 
as to possibly exhibit a series of quantized magnetization 
plateaux at m = 5(1), |(2), • • ■ , S + s — 1, its magne- 
tization curves have not systematically been studied so 
far. In spite of the vigorous argument, there are few re- 
ports on multi- plateau magnetization curves. It is true 
that a double-plateau structure lies in NH4CUCI3 pl| , p5[ , 
but it is owing to the variety of exchange interactions. 
We here demonstrate that the ferrimagnetic chain ( |1.2j ) 
generally exhibits a 2s -plateau magnetization curve with- 
out any anisotropy and any bond polymerization, namely, 
even at a = 1 and (5 = 0. We believe that the present cal- 
culations will accelerate physical measurements on vast 
ferrimagnetic chain compounds lying unexploited in the 
field of both inorganic and organic chemistry. 



The ground state of the isotropic Hamiltonian (1.2) 
without the Zeeman term, which is a multiplet of spin 
(S — s)N, exhibits elementary excitations of two distinct 
types jl(|. The excitations of ferromagnetic aspect, re- 
ducing the ground-state magnetization, form a gapless 
dispersion relation, whereas those of antifcrromagnetic 
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| p6[ . Considering that any magnetization plateau origi- 
nates from an antiferromagnetic interaction, it is convinc- 
ing that 2s of the total spin amplitude 5 + s contributes 
to the plateaux appearing. 

Let us turn back to Fig. |l] and observe the plateaux 
more carefully, especially as functions of a. In the 
cases of s = h, the plateaux are quite tough against 
the XY-\\\e anisotropy. They are stable all over the 
antifcrromagnetic-coupling region. These observations 
are in contrast with the classical behavior. The classi- 
cal spin- (S,s) Heisenberg Hamiltonian also exhibits the 
magnetization plateau at m = 5 — s, but it survives 
only a small amount of XY-like anisotropy. For in- 
stance, the critical value for (5, s) = (1, ^) is estimated as 
a c = 0.943(1) p| . The contrast between quantum spins 
and classical vectors suggests that the single plateaux 
in the quantum-spin magnetization curves may be at- 
tributed to the valence-bond excitation gap (valence-bond 
gap) rather than the N eel-state excitation gap (Neel gap). 
From this point of view, the a-dependences of the two 
coexistent plateaux in the cases of s = 1 are interesting. 
The tiny second plateaux are much more stable against 
the XY-like anisotropy than the steady-looking initial 
steps. Since the Neel state reaches the saturation via 
a single-step excitation, the second and higher plateaux 
should originate from the valence-bond gap. The magni- 
tude of the gap exponentially decreases with the increase 
of m, but the quantum gap-generation mechanism itself is 
rather tough against the XF-like anisotropy. The initial 
plateaux in the multi-step magnetization curves, which 
are relatively unstable against the XF-like anisotropy, 
may be attributed to the Neel gap. It seems that the 
lowest-lying magnetization plateaux are of quantum ap- 
pearance for s = i but are of classical aspect for s > 1. 
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FIG. 3. The variational ground-state phase diagrams on 
the SH-plane for the Hamiltonian (2) with a = 1 at various 
values of (S, s). 



In order to characterize the plateaux, we introduce 
a varia tion al wave function for the ground state of the 
model (pT3) as 



\g) = c N Yi(A]r s ( b }r\o) 

2s N 

+ E c vL IM) 2S ~V1) 2 ^(4 & ] - B My\V » 0-2) 



i=o j=i 



where cn and Cy B are the mixing coefficients. Since 
all the elemental states are asymptotically orthogonal to 
each other, the thermodynamic-limit variational ground 
states are considerably simple, as shown in Fig. ||, where 
we consider the Heisenberg point. The phase diagrams 
are exact on the line of S = 1, where the spin-(|, 1) chain, 
for example, reaches the saturation (S) via the double- 
bond dimer (DBD) and single-bond dimer (SBD) states. 
The point is that the 6 = ground states are better ap- 
proximated by the decoupled-dimer states than by the 
Neel (N) states in the cases of s — |, while vice versa in 
all o ther cases. For s — h, the variational wave function 
( |3.2|) ends up with cn = all over the S-H plain. The 
Neel-dimer crossover point S c is given by 



Sr. = 



2Ss - A + B + C 
A-B + C 



(3.3) 



with 



i=i 
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plateaux, with a quantum base, are stable over the whole 
antiferromagnetic-coupling region (a,b), while the initial 
plateaux in the multi-step process, taking on a classical 
character, less survive the XY-\ike anisotropy than the 
quantum higher plateaux (c,d). The existence of the sec- 
ond plateaux without any bond polymerization, which is 
the main issue in the present article, should be verified 
very carefully. So then we further employ an idea of the 
level spectroscopy thus called, in analyzing the sec- 
ond plateau. Comparing the relevant excitation energies 
whose scaling dimensions arc 2 at the critical point, we 
recognize the level crossing of them as the phase bound- 
ary. The thus-detected transitions are also plotted by 
broken lines in Fig. g. The slight difference between 
the two estimates inevitably a rise s from the logarithmic 
corrections to the scaling law (3.5), where the level spec- 
troscopy is more reliable than the naivest scaling analysis. 
Anyway, we may now fully be convinced of the existence 
of the novel multi-plateau magnetization curves. 



M' = Ni (s = 1), are encouraged. The plateaux would 
more or less be obscured in any actual measurement, but 
they, however small, should necessarily be detected by 
some anomaly in the magnetic susceptibility. The chem- 
ical modification of the bond alternation 5 and/or the 
exchange anisotropy a must help us to directly observe 
the second-step plateaux, though we take main interest 
in the Heisenberg point. 
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IV. SUMMARY AND FUTURE ASPECT 



The one-dimensional Heisenberg ferrimagnet with al- 
ternating spins S and s exhibits a 2s-plateau mag- 
netization curve even at the most symmetric point. 
It is interesting to compare the present spin-(S', s) 
ferrimagnetic chain with the spin-i bond-polymerized 
chain of period 2(S + s), that is, the (2S — 
l)-times-ferromagnetic-antiferromagnetic-(2s — l)-times- 
ferromagnetic-antiferromagnetic chain. In the strong 
ferromagnetic-coupling limit, the latter may be regarded 
as equivalent to the former. In the same meaning, 
the spin-iS antiferromagnetic Heisenberg chain can be 
viewed as a spin-^ bond-polymerized chain of period 
25. Such replica chains generally exhibit magnetiza- 
tion plateaux in certain regions of the ratio of the fer- 
romagnetic coupling Jp to the antiferromagnetic one 
Ja, 7 = Jf/Ja- However, in the case of the spin- 
i ferromagnetic-ferromagnetic-antiferromagnetic chain 
which is the replica model of the spin-| antiferromagnetic 
Heisenberg chain, Okamoto || and Hida |7]] reported that 
the plateau vanishes at 7 = 4 ~ 5 and therefore the pure 
Heisenberg chain exhibits no plateau. Thus the plateaux 
observed here in the most symmetric Heisenberg ferri- 
magnets still interest us to a great extent. 

As we come up the steps, the plateau length exponen- 
tially decreases. It is hard to numerically observe the 
higher-lying plateaux, still harder experimentally. Only 
the first and second plateaux may lie within the limits of 
measurement. In this context, we are fortunate to have 
a series of bimetallic quasi-one-dimensional complexes 
MM'(EDTA)-6H 2 (M, M' = Mn, Co, Ni, Cu) |§. Their 
exchange coupling constants are all about 10/cb[K] and 
thus the complete magnetization curves could techni- 
cally be observed. Magnetization measurements on them, 



especially with M = Mn(S 



J),Go(5 = §) and 
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TABLE I. Estimates of the antiferromagnetic excitation 
gap from the ground state for the Hamiltonian (2) with a — 1 
and 5 — by the use of the linear spin-wave theory (spin 
wave), the perturbation from the decoupled-dimer limit (per- 
turbation), and the exact diagonalization (exact). 



(S,s) 



spin wave 



perturbation 



exact 



(1.1) 

\2' 2> 
(1,1) 

(2,1) 
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1.759(1) 
2.842(1) 
1.615(10) 
2.730(5) 
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